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2.  Experimental  effort:  Visualization  of  oscillating  temperature  and  fiow  fields  (C.  Her¬ 
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Brief  description  of  projects 

Project  L  The  only  theory  presently  available  for  the  analysis  of  thermoacoustic  devices 
is  formulated  in  the  frequency  domain  and  is  only  applicable  in  the  linear  regime.  Its 
extension  to  the  nonlinear  case  appears  to  be  formidable  and  it  is  likely  that  nothing  short 
of  fully  numerical  simulations  will  be  available  for  a  quantitatively  accurate  application 
to  large-amplitude  pressure  oscillations.  For  this  reason,  an  approximate,  but  much  more 
manageable  model,  might  be  useful  both  as  an  exploratory  tool  and  as  a  design  guide.  The 
project  consists  in  the  development  of  such  a  model. 

A  small  amount  of  effort  has  also  been  spent  on  the  investigation  of  a  new  mechanism 
for  the  explanation  of  sonoluminescence. 

Project  2.  In  order  to  predict  heat  transfer  rates  in  the  heat  exchangers  of  thermoacoustic 
devices,  a  better  understanding  of  the  mechanisms  underlying  heat  transfer  in  oscillatory 
flows  is  essential.  We  investigate  such  physical  situations  using  holographic  interferometry 
(HI)  combined  with  high  speed  cinematography.  In  order  to  apply  HI  to  acoustically  driven 
flow  it  was  necessary  to  develop  a  new  evaluation  procedure  that  accounts  for  the  acoustic 
pressure  fluctuations.  This  evaluation  procedure  was  successfully  applied  to  measure  tem¬ 
perature  fields  as  well  as  heat  fluxes  in  a  thermoacoustic  refrigerator  model.  Current  efforts 
are  aiming  towards  the  application  of  HI  to  an  actual  thermoacoustic  refrigerator  that  has 
been  designed  applying  our  previously  developed  optimization  algorithm  based  on  the  short 
stack  boundary  layer  approximation. 

Project  3.  Analysis  and  prediction  of  the  behavior  of  thermoacoustic  devices  have  so  far 
relied  primarily  on  a  linearized  theory  based  on  simplified  quasi-one-dimensional  flow  mod¬ 
els.  In  order  to  overcome  these  restrictions,  this  effort  aims  at  (1)  the  development  of 
a  computational  methodology  for  the  simulation  of  the  non-linear  multi-dimensional  flow 
within  thermoacoustic  stacks  and  heat  exchangers,  and  (2)  the  implementation  of  the  re¬ 
sulting  numerical  schemes  to  analyze  the  fundamental  response  of  the  flow  within  critical 
components  and  to  quantify  the  performance  of  the  device.  The  computational  model  devel¬ 
opment  complements  parallel  efforts  aiming  at  gaining  detailed  fundamental  understanding 
of  thermoacoustic  devices  through  improved  analytical  models  and  laboratory  experiments. 
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ACCOMPLISHMENTS  -  PROJECT  1 


Time-domain  thermoacoustic  instability 


In  typical  thermoacoustic  prime  movers,  the  time  scale  for  the  development  of  the  insta¬ 
bility  is  long  compared  with  the  period  of  oscillation.  The  resulting  evolution  of  the  system 
acquires  the  character  of  modulated  oscillations  for  which  the  asymptotic  method  of  multiple 
time  scales  is  eminently  suitable. 

In  a  paper  recently  submitted  to  JASA  we  have  studied  this  problem  carrying  the  ex¬ 
pansion  to  second  order.  The  multiple-scales  results  can  be  further  simplified  by  using  the 
standard  so-called  short-stack  approximation.  We  have  found  that,  when  the  multiple-scales 
result  works  well,  so  does  the  short-stack  approximation.  The  final  expressions  have  therefore 
a  closed  analytical  form  that  enables  them  to  be  very  readily  evaluated. 

An  interesting  by-product  of  the  analysis  has  been  the  (approximate)  extension  of  the 
concept  of  critical  temperature  gradient  to  the  viscous  case.  The  result  is 


dTjr 


dx 


crit 


'up'  Ps  .  r-Tp  dPs/dx 
Cp  dPsjdx  ^  Ps 


Here  a  is  the  Prandtl  number,  uj  the  frequency,  Cp  the  specific  heat  at  constant  pressure, 
the  mean  temperature  in  the  stack,  and  P  =  P{x)  the  pressure  eigenfunction  of  the  mode 
of  interest;  the  subscript  S  indicates  evaluation  at  the  stack  position.  The  minimum  value 
of  this  critical  gradient  occurs  for 


dPs/dx 


(2) 


and  is 


dT,„ 


dx 


crit, min 


=  2  ^1  +  ^ 


Cp  ) 


1/2 


(3) 


In  the  paper  we  have  treated  not  only  the  case  of  a  tube  rigidly  terminated  at  both  ends, 
but  we  have  considered  all  combinations  of  closed/open  ends  and  all  modes  in  addition  to 
the  fundamental  one. 


Time-domain  approximation  of  the  exchange  terms 

It  will  be  recalled  from  our  previous  reports  that  the  model  developed  in  the  course  of 
this  work  has  been  obtained  by  integrating  the  conservation  equations  for  mass,  momentum, 
and  energy  over  the  cross  section  of  the  thermoacoustic  device  and  is,  accordingly,  quasi-one 
dimensional. 

In  last  year’s  report  we  pointed  out  that  the  necessity  to  use  approximate  expressions  for 
the  terms  describing  the  exchange  of  momentum  and  energy  between  the  stack  and  the  gas 
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worked  well  for  the  linear  theory,  but  gave  rise  to  some  problems  in  the  nonlinear  case.  In  our 
first  nonlinear  study  (submitted  to  JASA)  we  have  circumvented  the  problem  by  introducing 
ad  hoc  dissipation  terms  in  the  equations.  One  of  the  tasks  undertaken  since  then  has  been 
to  render  the  approach  more  systematic. 

The  exact  linear  relation  between  the  wall  shear  stress  r^u  and  the  mean  velocity  u  is: 


r  fv 

where  V  is  the  “wetted”  perimeter,  S  the  cross-sectional  area,  p  the  density,  and 

tanh  (1 -f 

"  {l+i)tl25v 


(4) 

(5) 


with  t  the  plate  spacing  and 

8v  = 

the  viscous  penetration  thickness  expressed  in  terms  of  the  kinematic  viscosity  v.  The 
relation  between  the  wall  heat  flux  and  the  mean  temperature  has  a  similar  structure. 

Our  nonlinear  model  requires  a  time-domain  form  for  the  frequency-domain  expression 
(4) .  In  principle,  such  a  form  can  be  derived  by  standard  Fourier  transform  theory,  but  takes 
on  the  form  of  a  convolution  integral  that  is  very  undesirable  for  numerical  purposes.  To 
tackle  the  problem  we  have  tried  to  approximate  (4)  by  a  form  of  the  type 

(1  -f-  iau))Tyj  =  p(b  +  icoj  —  u  (7) 

with  a,  b,  c,  d  suitably  chosen  real  constants.  In  the  time  domain,  this  expression  would 
correspond  to 

(l+a|)r.  =  /,(6  +  ic|  +  <i|j)«  (8) 

If  one  assumes  that,  in  a  low-order  sense,  u  satisfies  the  wave  equation  with  some  effective 
wave  speed  Cg,  an  alternative  form  of  (7)  would  be 

d 

dx^ 

We  wish  to  determine  the  constants  of  the  previous  expressions  not  by  a  numerical  fit, 
but  by  a  systematic  procedure.  We  have  tried  several  approaches,  and  in  the  end  we  have 
concluded  that  the  form  (9),  with  the  coefficients  determined  from  a  Taylor  series  expansion 
of  (4)  in  u)  around  the  frequency  of  the  fundamental  mode  gives  the  best  result.  We  are  in 
the  process  of  incorporating  this  new  approach  into  our  nonlinear  code  and  applying  it  to 
the  thermoacoustic  refrigerator  case. 


(1  -f  iauj)Tw  —  p  b  +  iao 
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Figure  1;  (a)  Time  dependent  temperature  measurements  with  the  conventional  approach 
and  the  required  corrections;  (b)  Contour  plot  of  the  time-averaged  temperature  distribution 
at  the  edge  of  the  thermoacoustic  refrigerator  model’s  stack  plates. 

ACCOMPLISHMENTS  -  PROJECT  2 


Temperature  and  heat  transfer  measurements  with  holographic  interferometry 
in  the  thermoacoustic  refrigerator  model. 

In  conventional  applications  of  holographic  interfereometry  (HI)  in  heat  transfer  mea¬ 
surements,  the  pressures  of  the  reference  state  and  of  the  measurement  state  are  maintained 
constant.  In  an  acoustically  driven  flow,  it  is  naturally  not  possible  to  maintain  the  mea¬ 
surement  state’s  pressure  constant.  In  order  to  resolve  this  problem  we  used  additional 
information  available  from  the  acoustic  field  and  developed  a  new  evaluation  procedure  that 
allows  temperature  measurements  with  HI  in  the  presence  of  pressure  variations.  However, 
we  have  found  that  the  conventional  approach  to  temperature  measurements  can  still  be 
applied  when  the  experimenter  is  interested  in  the  time  averaged  temperature  distribution 
(over  one  period  od  oscillation).  A  detailed  discussion  of  the  limitations  of  conventional 
temperature  measurements  with  HI  in  the  presence  of  pressure  variations  can  be  found  in 
the  recently  submitted  paper  by  Wetzel  and  Herman  (1997b). 

Figure  la  shows  the  time  dependent  temperature  measurements  at  a  fixed  point  within 
the  working  fluid  and  summarizes  some  of  the  important  aspects  of  HI  applied  to  an  acous¬ 
tically  driven  flow.  The  solid  circles  in  Figure  la  indicate  measurement  data  obtained  by  the 
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conventional  approach.  The  solid  line  indicates  a  correction  to  these  measurements  that  can 
be  calculated  from  the  error  function  derived  in  the  paper  by  Wetzel  and  herman  (1997b). 
Comparing  the  corrected  time-dependent  temperatures  with  the  measured  data,  indicated 
as  solid  squares  and  solid  circles  in  Figure  la,  respectively,  one  can  see  that  these  two  sets 
of  data  are  180  degrees  out  of  phase.  In  addition  to  the  phase  difference  we  c^n  also  see 
that  the  amplitude  of  the  small  temperature  fluctuation  decreases  from  about  4K  to  2K 
after  applying  the  described  correction.  When  neglecting  temperature  gradients  within  the 
working  fluid,  the  amplitude  of  the  small  temperature  fluctuations,  induced  by  the  acoustic 
oscillations,  can  be  roughly  estimated  by  a  simple  adiabatic  model  as 

(10) 

7  Pm 

In  this  equation,  Ta  denotes  the  amplitude  of  the  temperature  fluctuations,  7  the  ratio 
between  isobaric  and  isochoric  specific  heats,  Pa  the  peak  pressure  amplitude,  Pm  the  mean 
pressure  within  the  resonance  tube  and  the  mean  temperature.  Substituting  the  relevant 
values  for  our  experimental  conditions,  Tm  =  329  K,  Pa/Pm  ^  0.02,  and  7  =  1.4  into  Eq,  (10), 
we  obtain  an  amplitude  of  the  small  temperature  fluctuations  of  1.88  K.  This  value  agrees 
well  with  the  2K  of  the  corrected  temperature  distribution  shown  in  Fig.  la.  We  should 
also  note  that  both  the  measured  and  the  corrected  temperature  distributions  have  the  same 
value  for  the  mean  temperature,  as  indicated  by  the  dashed  line  in  Fig.  la.  From  this 
discussion,  we  can  conclude  that  the  conventional  approach  to  measure  temperatures  with 
HI  is  feasible  for  evaluation  of  time-averaged  temperature  distributions.  If  the  experimenter 
is  interested  in  the  time  resolved  temperature  distributions,  however,  it  is  necessary  to  add 
a  correction  as  described  by  Wetzel  and  Herman  (1997b). 

Figure  lb  shows  a  contour  plot  of  the  time-averaged  temperature  distribution  at  the  edge 
of  the  thermoacoustic  refrigerator  model’s  two  stack  plates  for  a  drive  ratio  of  2refrigerator 
model  is  also  described  in  the  paper  by  Wetzel  and  Herman  (1997b).  The  most  interesting 
feature  of  this  mean  temperature  distribution  is  the  shape  of  the  two  isotherms  that  enter  the 
lower  stack  plates  between  0  mm  <  a;  <  1  mm  and  are  marked  with  the  temperatures  57.50  °C 
and  58  ^C.  This  shape  indicates  that  the  mean  temperature  of  the  working  fluid  immediately 
above  the  stack  plates  is  higher  than  the  temperature  of  the  stack  plate.  Therefore  heat 
is  transferred  from  the  working  fluid  to  the  stack  plate.  This  behavior  is  different  when 
compared  to  a  steady  state  flow,  where  high  heat  transfer  rates  from  a  heated  surface  to 
the  working  fluid  are  achieved  at  the  entrance  of  a  parallel  plate  channel  because  of  the 
development  of  the  thermal  boundary  layer.  A  thorough  analysis  of  these  phenomena  is  one 
of  the  current  research  issues  and  results  will  be  published  in  the  near  future. 

Design  of  a  thermoacoustic  refrigerator  suitable  for  temperature  measurements 
with  HI 

Since  at  this  stage  of  our  research  the  problems  involved  in  the  application  of  HI  to 
an  acoustically  driven  flow  have  been  solved,  we  are  currently  working  on  measurements  of 
the  temperature  field  in  a  thermoacoustic  refrigerator.  The  design  of  this  thermoacoustic 
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Thermoacoustic 

refrigerator 

model 

Thermoacoustic 

refrigerator 

ThermoAcoustic 
Life  Sciences 
Refrigerator  TALSR 
(Garrett  1991) 

normalized  sthck 
length  ^  =  27rAx/A 

0.46 

0.14 

0.1188 

normalized  stack 
position  =  27rxc/ A 

0.78 

adjustable 

0.2 

normalized  plate 
spacing  5^h  = 

0.035 

0.15 

0.46 

blockage  ratio 

BR  =  h/ih  +  t) 

0.72 

0.83 

0.8 

Table  1:  Comparison  of  geometry  specific  design  parameters  between  the  thermoacoustic 
refrigerator  model  (Wetzel  and  Herman,  1997b),  the  thermoacoustic  refrigerator  (Brouillat, 
1996),  and  TALSR  Garrett  (1991). 

refrigerator  is  based  on  the  optimization  algorithm  developed  by  Wetzel  and  Herman  (1997a). 
One  advantage  of  this  algorithm  is  that  all  design  parameters  are  presented  in  normalized 
form.  Therefore,  it  allows  an  easy  comparison  between  different  devices. 

Table  1  shows  a  comparison  of  the  geometry  specific  design  parameters  between  the  ther¬ 
moacoustic  refrigerator  model  (Wetzel  and  Herman,  1997b),  the  new  thermoacoustic  refriger¬ 
ator  (Brouillat,  1996)  and  the  ThermoAcoustic  Life  Sciences  Refrigerator  (TALSR)  designed 
by  Garrett  (1991).  The  most  significant  improvement  when  comparing  the  thermoacoustic 
refrigerator  model  (previously  used  in  our  experiments)  with  the  new  thermoacoustic  refrig¬ 
erator  is  reflected  in  the  selection  of  the  plate  spacing.  In  the  current  experimental  setup  the 
plate  spacing  is  decreased  to  about  four  to  six  times  the  thermal  penetration  depth,  which 
results  in  a  normalized  plate  spacing  of  0.15.  Another  new  feature  of  the  thermoacoustic 
refrigerator  is  its  modular  structure  that  allows  easy  exchange  of  the  three  main  rnodules, 
(i)  stack,  (ii)  resonance  tube,  and  (iii)  loudspeaker.  This  structure  also  allows  an  adjustable 
stack  center  position.  When  comparing  the  geometry  specific  design  parameters  of  the  ther¬ 
moacoustic  refrigerator,  shown  in  Table  1,  with  the  values  used  in  the  design  of  TALSR,  one 
can  notice  the  good  agreement  between  the  two  designs.  A  detailed  description  of  the  new 
thermoacoustic  refrigerator’s  design  can  be  found  in  the  thesis  by  Brouillat  (1996). 

Figure  2  presents  a  snapshot  of  the  temperature  field  in  the  thermoacoustic  refrigerator. 
This  interferometric  image  was  recorded  with  the  high  speed  camera  at  a  framing  rate  of 
3000  frames  per  second.  The  black  rectangle  with  the  three  slots,  in  the  middle  of  Figure  2, 
corresponds  to  the  stack,  and  the  dark  and  bright  fringes  correspond  to  ’’quasi”-  isotherms, 
as  discussed  by  Wetzel  and  Herman  (1997b).  For  the  interferometric  image,  shown  in  Fig. 
2,  one  can  estimate  the  temperature  difference  from  a  dark  to  the  next  dark  fringe  to  be  of 
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Figure  2:  Snapshot  of  the  temperature  field  in  a  thermoacoustic  refrigerator. 

the  order  of  2K.  When  now  counting  the  dark  fringes  from  the  cold  side  to  the  hot  side,  we 
can  detect  5  dark  fringes  indicating  a  temperature  difference  of  about  lOK. 
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Figure  3:  Computed  values  of  the  temperature  difference  across  a  thermoacoustic  stack 
plotted  against  experimental  measurements  and  predictions  of  the  linear  theory  (adapted 
from  [8]). 


ACCOMPLISHMENTS  -  PROJECT  3 

Thermoacoustic  devices  are  typically  characterized  by  a  multitude  of  disparate  length- 
scales,  which  renders  analysis  or  simulation  of  the  non-linear  equations  of  motion  extremely 
difficult.  In  order  to  overcome  these  scale-related  difficulties,  we  have  been  involved  in  the 
development  of  an  efficient,  multi-dimensional  flow  solver  that  is  suitable  for  the  study  of 
thermo- fluid  phenomena  in  the  neighborhood  of  thermoacoustic  stacks.  The  early  version 
of  this  model  is  described  in  a  paper  recently  published  in  the  Jouvnal  of  CoTfiputatioTicil 
Physics  [1]. 

In  the  past,  we  have  focused  our  attention  on  stacks  operating  at  low  drive  ratios,  and  thus 
implemented  a  simplified  version  of  the  model  which  ignores  thermal  stratification  effects. 
The  simplified  model  has  enabled  us  to  analyze  the  essential  features  of  the  unsteady  flow, 
and  to  perform  a  detailed  analysis  of  the  behavior  of  hydrodynamic  (mechanical  energy) 
losses  within  the  stack  [2]. 

In  the  past  year,  we  have  generalized  the  original  model  in  order  to  account  for  thermal 
stratification  effects,  and  thus  extended  the  simulations  to  moderate/high  drive  ratios.  The 
extension  required  careful  design  and  analysis  of  a  new  fast  solver,  and  its  optimization  in 
order  to  maintain  the  efficiency  of  the  computations.  A  paper  describing  the  extended  model 
has  been  recently  submitted  for  publication  [3]. 

In  addition  to  this  model  development  effort,  we  have  performed  extensive  tests  aiming,  in 
particular,  at  validation  of  the  extended  computations  [4,  5].  To  this  end,  the  model  was  used 
to  compute  the  temperature  differences  across  a  thermoacoustic  stack  at  no  load  conditions, 


8 


and  the  results  were  compared  to  experimental  observations  of  Atchley  et  al.  [6].  A  sample 
of  the  computations  is  given  in  figure  3,  which  compares  the  computed  predictions  of  the 
temperature  drop  across  the  stack  to  the  experimental  observations  [6]  and  the  predictions 
of  the  linear  theory  [7]. 

With  the  extension  of  the  computations  to  stratified  flow,  and  their  validation  against 
experimental  data,  we  have  reached  an  important  milestone  in  the  development  of  the  multi¬ 
dimensional  flow  code.  Briefly,  we  are  now  at  stage  where  the  performance  of  simplified 
stacks  or  simplified  stacks/heat  exchanger  configurations  (figure  4  [8])  can  be  predicted  and 
analyzed  based  on  first  principles,  with  minimal  or  no  need  for  ad-hoc  modeling  assumptions. 


Cooling  load  (Qc) 

Figure  4;  Variation  of  relative  COP  with  the  cooling  load.  The  markers  represent  the  results 
of  multi-dimensional  computations  (adapted  from  [8]). 

In  the  coming  year,  we  plan  to  further  extend  the  current  stack  model  into  a  model  for 
the  entire  thermoacoustic  refrigerator.  To  this  end,  we  plan  to  couple  the  current  low-Mach- 
number  solver  with  a  quasi-lD  solver  for  weakly-nonlinear  resonance  tube  acoustics.  We 
plan  to  adapt  some  of  our  asymptotics-based  schemes  [9,  10]  for  this  purpose. 
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